Circuit-wide gene network analysis reveals sex-specific roles for phosphodiesterase 1b in cocaine addiction.

Cocaine use disorder is a significant public health issue without an effective pharmacological treatment. Successful treatments are hindered in part by an incomplete understanding of the molecular mechanisms that underlie long-lasting maladaptive plasticity and addiction-like behaviors. Here, we leverage a large RNA-sequencing dataset to generate gene co-expression networks across 6 interconnected regions of the brain's reward circuitry from mice that underwent saline or cocaine self-administration. We identify phosphodiesterase 1b (Pde1b), a Ca2+/calmodulin-dependent enzyme that increases cAMP and cGMP hydrolysis, as a central hub gene within a nucleus accumbens (NAc) gene module that was bioinformatically associated with addiction-like behavior. Chronic cocaine exposure increases Pde1b expression in NAc D2 medium spiny neurons (MSNs) in male but not female mice. Viral-mediated Pde1b overexpression in NAc reduces cocaine self-administration in female rats, but increases seeking in both sexes. In female mice, overexpressing Pde1b in D1 MSNs attenuates the locomotor response to cocaine, with the opposite effect in D2 MSNs. Overexpressing Pde1b in D1/D2 MSNs had no effect on the locomotor response to cocaine in male mice. At the electrophysiological level, Pde1b overexpression reduces sEPSC frequency in D1 MSNs, while increasing excitability of D2 MSNs. Lastly, Pde1b overexpression significantly reduced the number of differentially expressed genes (DEGs) in NAc following chronic cocaine, with discordant effects on gene transcription between sexes. Together, we identify novel gene modules across the brain's reward circuitry associated with addiction-like behavior and explore the role of Pde1b in regulating the molecular, cellular, and behavioral responses to cocaine.Significance Statement Cocaine use disorder is a major public health challenge without an effective pharmacological treatment. Here, we leverage a combination of genome-wide RNA sequencing, gene co-expression network analysis, and bioinformatic analyses of cocaine self-administration behavior to identify a role for phosphodiesterase 1b (Pde1b) in regulating maladaptive, addiction-like behavior. Our studies reveal cell-type- and sex-specific roles for Pde1b in regulating the molecular, cellular, and behavioral responses to cocaine, yielding insight into the molecular mechanisms by which cocaine induces maladaptive plasticity in the brain's reward circuity to drive addiction-like behavior. These discoveries guide directions for future research investigating the molecular basis of cocaine action and provide a pathway for therapeutic development for cocaine use disorder.

Transcriptomic landscape of human induced pluripotent stem cell-derived osteogenic differentiation identifies a regulatory role of KLF16.

Osteogenic differentiation is essential for bone development and metabolism, but the underlying gene regulatory networks have not been well investigated. We differentiated mesenchymal stem cells, derived from 20 human induced pluripotent stem cell lines, into preosteoblasts and osteoblasts, and performed systematic RNA-seq analyses of 60 samples for differential gene expression. We noted a highly significant correlation in expression patterns and genomic proximity among transcription factor (TF) and long noncoding RNA (lncRNA) genes. We identified TF-TF regulatory networks, regulatory roles of lncRNAs on their neighboring coding genes for TFs and splicing factors, and differential splicing of TF, lncRNA, and splicing factor genes. TF-TF regulatory and gene co-expression network analyses suggested an inhibitory role of TF KLF16 in osteogenic differentiation. We demonstrate that in vitro overexpression of human KLF16 inhibits osteogenic differentiation and mineralization, and in vivo Klf16+/- mice exhibit increased bone mineral density, trabecular number, and cortical bone area. Thus, our model system highlights the regulatory complexity of osteogenic differentiation and identifies novel osteogenic genes.

Multiscale protein networks systematically identify aberrant protein interactions and oncogenic regulators in seven cancer types.

Global proteomic data generated by advanced mass spectrometry (MS) technologies can help bridge the gap between genome/transcriptome and functions and hold great potential in elucidating unbiased functional models of pro-tumorigenic pathways. To this end, we collected the high-throughput, whole-genome MS data and conducted integrative proteomic network analyses of 687 cases across 7 cancer types including breast carcinoma (115 tumor samples; 10,438 genes), clear cell renal carcinoma (100 tumor samples; 9,910 genes), colorectal cancer (91 tumor samples; 7,362 genes), hepatocellular carcinoma (101 tumor samples; 6,478 genes), lung adenocarcinoma (104 tumor samples; 10,967 genes), stomach adenocarcinoma (80 tumor samples; 9,268 genes), and uterine corpus endometrial carcinoma UCEC (96 tumor samples; 10,768 genes). Through the protein co-expression network analysis, we identified co-expressed protein modules enriched for differentially expressed proteins in tumor as disease-associated pathways. Comparison with the respective transcriptome network models revealed proteome-specific cancer subnetworks associated with heme metabolism, DNA repair, spliceosome, oxidative phosphorylation and several oncogenic signaling pathways. Cross-cancer comparison identified highly preserved protein modules showing robust pan-cancer interactions and identified endoplasmic reticulum-associated degradation (ERAD) and N-acetyltransferase activity as the central functional axes. We further utilized these network models to predict pan-cancer protein regulators of disease-associated pathways. The top predicted pan-cancer regulators including RSL1D1, DDX21 and SMC2, were experimentally validated in lung, colon, breast cancer and fetal kidney cells. In summary, this study has developed interpretable network models of cancer proteomes, showcasing their potential in unveiling novel oncogenic regulators, elucidating underlying mechanisms, and identifying new therapeutic targets.

Shared and divergent transcriptomic regulation in nucleus accumbens D1 and D2 medium spiny neurons by cocaine and morphine.

Substance use disorders (SUDs) induce widespread molecular dysregulation in the nucleus accumbens (NAc), a brain region pivotal for coordinating motivation and reward. These molecular changes are thought to support lasting neural and behavioral disturbances that promote drug-seeking in addiction. However, different drug classes exert unique influences on neural circuits, cell types, physiology, and gene expression despite the overlapping symptomatology of SUDs. To better understand common and divergent molecular mechanisms governing SUD pathology, our goal was to survey cell-type-specific restructuring of the NAc transcriptional landscape in after psychostimulant or opioid exposure. We combined fluorescence-activated nuclei sorting and RNA sequencing to profile NAc D1 and D2 medium spiny neurons (MSNs) across cocaine and morphine exposure paradigms, including initial exposure, prolonged withdrawal after repeated exposure, and re-exposure post-withdrawal. Our analyses reveal that D1 MSNs display many convergent transcriptional responses across drug classes during exposure, whereas D2 MSNs manifest mostly divergent responses between cocaine and morphine, with morphine causing more adaptations in this cell type. Utilizing multiscale embedded gene co-expression network analysis (MEGENA), we discerned transcriptional regulatory networks subserving biological functions shared between cocaine and morphine. We observed largely integrative engagement of overlapping gene networks across drug classes in D1 MSNs, but opposite regulation of key D2 networks, highlighting potential therapeutic gene network targets within MSNs. These studies establish a landmark, cell-type-specific atlas of transcriptional regulation induced by cocaine and by morphine that can serve as a foundation for future studies towards mechanistic understanding of SUDs. Our findings, and future work leveraging this dataset, will pave the way for the development of targeted therapeutic interventions, addressing the urgent need for more effective treatments for cocaine use disorder and enhancing the existing strategies for opioid use disorder.

Smarca4 deficiency induces Pttg1 oncogene upregulation and hyperproliferation of tubular and interstitial cells during kidney development.

Kidney formation and nephrogenesis are controlled by precise spatiotemporal gene expression programs, which are coordinately regulated by cell-cycle, cell type-specific transcription factors and epigenetic/chromatin regulators. However, the roles of epigenetic/chromatin regulators in kidney development and disease remain poorly understood. In this study, we investigated the impact of deleting the chromatin remodeling factor Smarca4 (Brg1), a human Wilms tumor-associated gene, in Wnt4-expressing cells. Smarca4 deficiency led to severe tubular defects and a shortened medulla. Through unbiased single-cell RNA sequencing analyses, we identified multiple types of Wnt4 Cre-labeled interstitial cells, along with nephron-related cells. Smarca4 deficiency increased interstitial cells but markedly reduced tubular cells, resulting in cells with mixed identity and elevated expression of cell-cycle regulators and genes associated with extracellular matrix and epithelial-to-mesenchymal transition/fibrosis. We found that Smarca4 loss induced a significant upregulation of the oncogene Pttg1 and hyperproliferation of Wnt4 Cre-labeled cells. These changes in the cellular state could hinder the cellular transition into characteristic tubular structures, eventually leading to fibrosis. In conclusion, our findings shed light on novel cell types and genes associated with Wnt4 Cre-labeled cells and highlight the critical role of Smarca4 in regulating tubular cell differentiation and the expression of the cancer-causing gene Pttg1 in the kidney. These findings may provide valuable insights into potential therapeutic strategies for renal cell carcinoma resulting from SMARCA4 deficiency.

Molecular differences in brain regional vulnerability to aging between males and females.

Background: Aging-related cognitive decline is associated with brain structural changes and synaptic loss. However, the molecular mechanisms of cognitive decline during normal aging remain elusive. Results: Using the GTEx transcriptomic data from 13 brain regions, we identified aging-associated molecular alterations and cell-type compositions in males and females. We further constructed gene co-expression networks and identified aging-associated modules and key regulators shared by both sexes or specific to males or females. A few brain regions such as the hippocampus and the hypothalamus show specific vulnerability in males, while the cerebellar hemisphere and the anterior cingulate cortex regions manifest greater vulnerability in females than in males. Immune response genes are positively correlated with age, whereas those involved in neurogenesis are negatively correlated with age. Aging-associated genes identified in the hippocampus and the frontal cortex are significantly enriched for gene signatures implicated in Alzheimer's disease (AD) pathogenesis. In the hippocampus, a male-specific co-expression module is driven by key synaptic signaling regulators including VSNL1, INA, CHN1 and KCNH1; while in the cortex, a female-specific module is associated with neuron projection morphogenesis, which is driven by key regulators including SRPK2, REPS2 and FXYD1. In the cerebellar hemisphere, a myelination-associated module shared by males and females is driven by key regulators such as MOG, ENPP2, MYRF, ANLN, MAG and PLP1, which have been implicated in the development of AD and other neurodegenerative diseases. Conclusions: This integrative network biology study systematically identifies molecular signatures and networks underlying brain regional vulnerability to aging in males and females. The findings pave the way for understanding the molecular mechanisms of gender differences in developing neurodegenerative diseases such as AD.

Sex specific molecular networks and key drivers of Alzheimer's disease.

BACKGROUND: Alzheimer's disease (AD) is a progressive and age-associated neurodegenerative disorder that affects women disproportionally. However, the underlying mechanisms are poorly characterized. Moreover, while the interplay between sex and ApoE genotype in AD has been investigated, multi-omics studies to understand this interaction are limited. Therefore, we applied systems biology approaches to investigate sex-specific molecular networks of AD. METHODS: We integrated large-scale human postmortem brain transcriptomic data of AD from two cohorts (MSBB and ROSMAP) via multiscale network analysis and identified key drivers with sexually dimorphic expression patterns and/or different responses to APOE genotypes between sexes. The expression patterns and functional relevance of the top sex-specific network driver of AD were further investigated using postmortem human brain samples and gene perturbation experiments in AD mouse models. RESULTS: Gene expression changes in AD versus control were identified for each sex. Gene co-expression networks were constructed for each sex to identify AD-associated co-expressed gene modules shared by males and females or specific to each sex. Key network regulators were further identified as potential drivers of sex differences in AD development. LRP10 was identified as a top driver of the sex differences in AD pathogenesis and manifestation. Changes of LRP10 expression at the mRNA and protein levels were further validated in human AD brain samples. Gene perturbation experiments in EFAD mouse models demonstrated that LRP10 differentially affected cognitive function and AD pathology in sex- and APOE genotype-specific manners. A comprehensive mapping of brain cells in LRP10 over-expressed (OE) female E4FAD mice suggested neurons and microglia as the most affected cell populations. The female-specific targets of LRP10 identified from the single cell RNA-sequencing (scRNA-seq) data of the LRP10 OE E4FAD mouse brains were significantly enriched in the LRP10-centered subnetworks in female AD subjects, validating LRP10 as a key network regulator of AD in females. Eight LRP10 binding partners were identified by the yeast two-hybrid system screening, and LRP10 over-expression reduced the association of LRP10 with one binding partner CD34. CONCLUSIONS: These findings provide insights into key mechanisms mediating sex differences in AD pathogenesis and will facilitate the development of sex- and APOE genotype-specific therapies for AD.

Multiscale network analysis identifies potential receptors for SARS-CoV-2 and reveals their tissue-specific and age-dependent expression.

The coronavirus disease 2019 (COVID-19) pandemic has affected tens of millions of individuals and caused hundreds of thousands of deaths worldwide. Here, we present a comprehensive, multiscale network analysis of the transcriptional response to the virus. In particular, we focused on key regulators, cell receptors, and host processes that were hijacked by the virus for its advantage. ACE2-controlled processes involved CD300e (a TYROBP receptor) as a key regulator and the activation of IL-2 pro-inflammatory cytokine signaling. We further investigated the age dependency of such receptors in different tissues. In summary, this study provides novel insights into the gene regulatory organization during the SARS-CoV-2 infection and the tissue-specific, age-dependent expression of the cell receptors involved in COVID-19.

Pseudo-temporal dynamics of chemoresistant triple negative breast cancer cells reveal EGFR/HER2 inhibition as synthetic lethal during mid-neoadjuvant chemotherapy.

In the absence of targetable hormonal axes, chemoresistance for triple-negative breast cancer (TNBC) often compromises patient outcomes. To investigate the underlying tumor dynamics, we performed trajectory analysis on the single-nuclei RNA-seq (snRNA-seq) of chemoresistant tumor clones during neoadjuvant chemotherapy (NAC). It revealed a common tumor trajectory across multiple patients with HER2-like expansions during NAC. Genome-wide CRISPR-Cas9 knock-out on mammary epithelial cells revealed chemosensitivity-promoting knock-outs were up-regulated along the tumor trajectory. Furthermore, we derived a consensus gene signature of TNBC chemoresistance by comparing the trajectory transcriptome with chemoresistant transcriptomes from TNBC cell lines and poor prognosis patient samples to predict FDA-approved drugs, including afatinib (pan-HER inhibitor), targeting the consensus signature. We validated the synergistic efficacy of afatinib and paclitaxel in chemoresistant TNBC cells and confirmed pharmacological suppression of the consensus signature. The study provides a dynamic model of chemoresistant tumor transcriptome, and computational framework for pharmacological intervention.

Multi-omics-based analysis of high grade serous ovarian cancer subtypes reveals distinct molecular processes linked to patient prognosis.

Despite advancements in treatment, high-grade serous ovarian cancer (HGSOC) is still characterized by poor patient outcomes. To understand the molecular heterogeneity of this disease, which underlies the challenge in selecting optimal treatments for HGSOC patients, we have integrated genomic, transcriptomic, and epigenetic information to identify seven new HGSOC subtypes using a multiscale clustering method. These subtypes not only have significantly distinct overall survival, but also exhibit unique patterns of gene expression, microRNA expression, DNA methylation, and copy number alterations. As determined by our analysis, patients with similar clinical outcomes have distinct profiles of activated or repressed cellular processes, including cell cycle, epithelial-to-mesenchymal transition, immune activation, interferon response, and cilium organization. Furthermore, we performed a multiscale gene co-expression network analysis to identify subtype-specific key regulators and predicted optimal targeted therapies based on subtype-specific gene expression. In summary, this study provides new insights into the cellular heterogeneity of the HGSOC genomic, epigenetic, and transcriptomic landscapes and provides a basis for future studies into precision medicine for HGSOC patients.

Common and species-specific molecular signatures, networks, and regulators of influenza virus infection in mice, ferrets, and humans.

Molecular responses to influenza A virus (IAV) infections vary between mammalian species. To identify conserved and species-specific molecular responses, we perform a comparative study of transcriptomic data derived from blood cells, primary epithelial cells, and lung tissues collected from IAV-infected humans, ferrets, and mice. The molecular responses in the human host have unique functions such as antigen processing that are not observed in mice or ferrets. Highly conserved gene coexpression modules across the three species are enriched for IAV infection-induced pathways including cell cycle and interferon (IFN) signaling. TDRD7 is predicted as an IFN-inducible host factor that is up-regulated upon IAV infection in the three species. TDRD7 is required for antiviral IFN response, potentially modulating IFN signaling via the JAK/STAT/IRF9 pathway. Identification of the common and species-specific molecular signatures, networks, and regulators of IAV infection provides insights into host-defense mechanisms and will facilitate the development of novel therapeutic interventions against IAV infection.

Crystallin Mu in Medial Amygdala Mediates the Effect of Social Experience on Cocaine Seeking in Males but Not in Females.

BACKGROUND: Social experiences influence susceptibility to substance use disorder. The adolescent period is associated with the development of social reward and is exceptionally sensitive to disruptions to reward-associated behaviors by social experiences. Social isolation (SI) during adolescence alters anxiety- and reward-related behaviors in adult males, but little is known about females. The medial amygdala (meA) is a likely candidate for the modulation of social influence on drug reward because it regulates social reward, develops during adolescence, and is sensitive to social stress. However, little is known regarding how the meA responds to drugs of abuse. METHODS: We used adolescent SI coupled with RNA sequencing to better understand the molecular mechanisms underlying meA regulation of social influence on reward. RESULTS: We show that SI in adolescence, a well-established preclinical model for addiction susceptibility, enhances preference for cocaine in male but not in female mice and alters cocaine-induced protein and transcriptional profiles within the adult meA particularly in males. To determine whether transcriptional mechanisms within the meA are important for these behavioral effects, we manipulated Crym expression, a sex-specific key driver gene identified through differential gene expression and coexpression network analyses, specifically in meA neurons. Overexpression of Crym, but not another key driver that did not meet our sex-specific criteria, recapitulated the behavioral and transcriptional effects of adolescent SI. CONCLUSIONS: These results show that the meA is essential for modulating the sex-specific effects of social experience on drug reward and establish Crym as a critical mediator of sex-specific behavioral and transcriptional plasticity.

Dynamic plasticity of prostate cancer intermediate cells during androgen receptor-targeted therapy.

Treatment-emergent small cell neuroendocrine prostate cancer (t-SCNC) is associated with an epithelial lineage switch from an androgen receptor (AR)-positive to neuroendocrine (NE)-marker-positive status. Understanding the potential for reversibility of this aggressive disease state has been hampered by the paucity of models suitable for studying rate-limiting, transitional, or intermediate tumor cell subpopulations. We define a dual reporter model that measures acute transcriptional changes in response to castration or AR targeting agents. We identify steady-state transcriptional heterogeneity in AR and NE biomarkers, including intermediate subpopulations that are coordinately high for prostate-specific antigen (PSA) and neuron-specific enoclase (NSE) promoter activity. In the presence of castration or AR inhibitors, intermediate cells were necessary and sufficient for therapy-induced conversion of human PC cells to an NSE-high transcriptional status. Using hormone add-back studies, treatment-induced PSA-NSE transcriptional plasticity was reversible in PTEN-deficient PC cells but not in the presence of secondary genetic driver genes, including MYCN.

Guidelines for bioinformatics of single-cell sequencing data analysis in Alzheimer's disease: review, recommendation, implementation and application.

Alzheimer's disease (AD) is the most common form of dementia, characterized by progressive cognitive impairment and neurodegeneration. Extensive clinical and genomic studies have revealed biomarkers, risk factors, pathways, and targets of AD in the past decade. However, the exact molecular basis of AD development and progression remains elusive. The emerging single-cell sequencing technology can potentially provide cell-level insights into the disease. Here we systematically review the state-of-the-art bioinformatics approaches to analyze single-cell sequencing data and their applications to AD in 14 major directions, including 1) quality control and normalization, 2) dimension reduction and feature extraction, 3) cell clustering analysis, 4) cell type inference and annotation, 5) differential expression, 6) trajectory inference, 7) copy number variation analysis, 8) integration of single-cell multi-omics, 9) epigenomic analysis, 10) gene network inference, 11) prioritization of cell subpopulations, 12) integrative analysis of human and mouse sc-RNA-seq data, 13) spatial transcriptomics, and 14) comparison of single cell AD mouse model studies and single cell human AD studies. We also address challenges in using human postmortem and mouse tissues and outline future developments in single cell sequencing data analysis. Importantly, we have implemented our recommended workflow for each major analytic direction and applied them to a large single nucleus RNA-sequencing (snRNA-seq) dataset in AD. Key analytic results are reported while the scripts and the data are shared with the research community through  GitHub. In summary, this comprehensive review provides insights into various approaches to analyze single cell sequencing data and offers specific guidelines for study design and a variety of analytic directions. The review and the accompanied software tools will serve as a valuable resource for studying cellular and molecular mechanisms of AD, other diseases, or biological systems at the single cell level.

Lifelong chronic psychosocial stress induces a proteomic signature of Alzheimer's disease in wildtype mice.

Late onset, sporadic Alzheimer's disease (AD) accounts for the vast majority of cases. Unlike familial AD, the factors that drive the onset of sporadic AD are poorly understood, although aging and stress play a role. The early onset/severity of neuropathology observed in most genetic mouse models of AD hampers the study of the role of aging and environmental factors; thus alternate strategies are necessary to understand the contributions of these factors to sporadic AD. We demonstrate that mice acquiring a low social status (subordinate) in a lifelong chronic psychosocial stress (CPS) model, accrue widespread proteomic changes in the frontal/temporal cortex during aging. To better understand the significance of these stress-induced changes, we compared the differentially expressed proteins (DEPs) of subordinate mice to those of patients at varying stages of dementia. Sixteen and fifteen DEPs upregulated in subordinate mice were also upregulated in patients with mild cognitive impairment (MCI) and AD, respectively. Six of those upregulated proteins (CPE, ERC2, GRIN2B, SLC6A1, SYN1, WFS1) were shared by subordinate mice and patients with MCI or AD. Finally, comparison with a spatially detailed transcriptomic database revealed that the superior frontal gyrus and hippocampus had the greatest overlap between mice subjected to lifelong CPS and AD patients. Overall, most of the overlapping proteins were functionally associated with enhanced NMDA receptor mediated glutamatergic signaling, an excitotoxicity mechanism known to affect neurodegeneration. These findings support the association between stress and AD progression and provide valuable insight into potential early biomarkers and protein mediators of this relationship.

Sex-Specific Transcriptional Changes in Response to Adolescent Social Stress in the Brain's Reward Circuitry.

BACKGROUND: Sex differences in addiction have been described in humans and animal models. A key factor that influences addiction in both males and females is adolescent experience. Adolescence is associated with higher vulnerability to substance use disorders, and male rodents subjected to adolescent social isolation (SI) stress form stronger preferences for drugs of abuse in adulthood. However, little is known about how females respond to SI, and few studies have investigated the transcriptional changes induced by SI in the brain's reward circuitry. METHODS: We tested the hypothesis that SI alters the transcriptome in a persistent and sex-specific manner in prefrontal cortex, nucleus accumbens, and ventral tegmental area. Mice were isolated or group housed from postnatal day P22 to P42, then group housed until ∼P90. Transcriptome-wide changes were investigated by RNA sequencing after acute or chronic cocaine or saline administration. RESULTS: We found that SI disrupts sex-specific transcriptional responses to cocaine and reduces sex differences in gene expression across all three brain regions. Furthermore, SI induces gene expression profiles in males that more closely resemble group-housed females, suggesting that SI "feminizes" the male transcriptome. Coexpression analysis reveals that such disruption of sex differences in gene expression alters sex-specific gene networks and identifies potential sex-specific key drivers of these transcriptional changes. CONCLUSIONS: Together, these data show that SI has region-specific effects on sex-specific transcriptional responses to cocaine and provide a better understanding of reward-associated transcription that differs in males and females.

Single-cell analysis identifies a key role for Hhip in murine coronal suture development.

Craniofacial development depends on formation and maintenance of sutures between bones of the skull. In sutures, growth occurs at osteogenic fronts along the edge of each bone, and suture mesenchyme separates adjacent bones. Here, we perform single-cell RNA-seq analysis of the embryonic, wild type murine coronal suture to define its population structure. Seven populations at E16.5 and nine at E18.5 comprise the suture mesenchyme, osteogenic cells, and associated populations. Expression of Hhip, an inhibitor of hedgehog signaling, marks a mesenchymal population distinct from those of other neurocranial sutures. Tracing of the neonatal Hhip-expressing population shows that descendant cells persist in the coronal suture and contribute to calvarial bone growth. In Hhip-/- coronal sutures at E18.5, the osteogenic fronts are closely apposed and the suture mesenchyme is depleted with increased hedgehog signaling compared to those of the wild type. Collectively, these data demonstrate that Hhip is required for normal coronal suture development.

Prosaposin mediates inflammation in atherosclerosis.

Macrophages play a central role in the pathogenesis of atherosclerosis. The inflammatory properties of these cells are dictated by their metabolism, of which the mechanistic target of rapamycin (mTOR) signaling pathway is a key regulator. Using myeloid cell-specific nanobiologics in apolipoprotein E-deficient (Apoe -/-) mice, we found that targeting the mTOR and ribosomal protein S6 kinase-1 (S6K1) signaling pathways rapidly diminished plaque macrophages' inflammatory activity. By investigating transcriptome modifications, we identified Psap, a gene encoding the lysosomal protein prosaposin, as closely related with mTOR signaling. Subsequent in vitro experiments revealed that Psap inhibition suppressed both glycolysis and oxidative phosphorylation. Transplantation of Psap -/- bone marrow to low-density lipoprotein receptor knockout (Ldlr -/-) mice led to a reduction in atherosclerosis development and plaque inflammation. Last, we confirmed the relationship between PSAP expression and inflammation in human carotid atherosclerotic plaques. Our findings provide mechanistic insights into the development of atherosclerosis and identify prosaposin as a potential therapeutic target.

Diversified transcriptional responses of myeloid and glial cells in spinal cord injury shaped by HDAC3 activity.

The innate immune response influences neural repair after spinal cord injury (SCI). Here, we combined myeloid-specific transcriptomics and single-cell RNA sequencing to uncover not only a common core but also temporally distinct gene programs in injury-activated microglia and macrophages (IAM). Intriguingly, we detected a wide range of microglial cell states even in healthy spinal cord. Upon injury, IAM progressively acquired overall reparative, yet diversified transcriptional profiles, each comprising four transcriptional subtypes with specialized tasks. Notably, IAM have both distinct and common gene signatures as compared to neurodegeneration-associated microglia, both engaging phagocytosis, autophagy, and TyroBP pathways. We also identified an immediate response microglia subtype serving as a source population for microglial transformation and a proliferative subtype controlled by the epigenetic regulator histone deacetylase 3 (HDAC3). Together, our data unveil diversification of myeloid and glial subtypes in SCI and an extensive influence of HDAC3, which may be exploited to enhance functional recovery.

Network models of primary melanoma microenvironments identify key melanoma regulators underlying prognosis.

Melanoma is the most lethal skin malignancy, driven by genetic and epigenetic alterations in the complex tumour microenvironment. While large-scale molecular profiling of melanoma has identified molecular signatures associated with melanoma progression, comprehensive systems-level modeling remains elusive. This study builds up predictive gene network models of molecular alterations in primary melanoma by integrating large-scale bulk-based multi-omic and single-cell transcriptomic data. Incorporating clinical, epigenetic, and proteomic data into these networks reveals key subnetworks, cell types, and regulators underlying melanoma progression. Tumors with high immune infiltrates are found to be associated with good prognosis, presumably due to induced CD8+ T-cell cytotoxicity, via MYO1F-mediated M1-polarization of macrophages. Seventeen key drivers of the gene subnetworks associated with poor prognosis, including the transcription factor ZNF180, are tested for their pro-tumorigenic effects in vitro. The anti-tumor effect of silencing ZNF180 is further validated using in vivo xenografts. Experimentally validated targets of ZNF180 are enriched in the ZNF180 centered network and the known pathways such as melanoma cell maintenance and immune cell infiltration. The transcriptional networks and their critical regulators provide insights into the molecular mechanisms of melanomagenesis and pave the way for developing therapeutic strategies for melanoma.